Ferromagnetic resonance force spectroscopy of individual sub-micron size samples 
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We review how a magnetic resonance force microscope (MRFM) can be applied to perform fer- 
romagnetic resonance (FMR) spectroscopy of individual sub-micron size samples. We restrict our 
attention to a thorough study of the spin-wave eigen-modes excited in permalloy (Py) disks pat- 
terned out of the same 43.3 nm thin film. The disks have a diameter of either 1.0 or 0.5 /im and 
are quasi-saturated by a perpendicularly applied magnetic field. It is shown that quantitative spec- 
troscopic information can be extracted from the MRFM measurements. In particular, the data are 
extensively compared with complementary approximate models of the dynamical susceptibility: i) 
a 2D analytical model, which assumes an homogeneous magnetization dynamics along the thickness 
and ii) a full 3D micromagnetic simulation, which assumes an homogeneous magnetization dynamics 
below a characteristic length scale c and which approximates the cylindrical sample volume by a 
discretized representation with regular cubic mesh of lateral size c = 3.9 nm. In our analysis, the dis- 
tortions due to a breaking of the axial symmetry are taken into account, both models incorporating 
the possibility of a small misalignment between the applied field and the normal of the disks. 



I. INTRODUCTION 

Development of innovative tools capable of measur- 
ing the local magnetization dynamics M(t,r) inside a 
ferromagnetic nano-structure is an objective of primary 
importance. New technology fields related to magnetic 
materials, like spintronics, depend on one's ability to an- 
alyze and predict the out-of-equilibrium state of M in 
nanoscalc hybrid structures. In this respect, models de- 
rived from microscopic principles are often limited be- 
cause of the numerous degrees of freedom strongly cou- 
pled to the magnetization. Hence, phenomcnological ap- 
proaches are often used, and require a constant compar- 
ison with experiments. 

Several original techniques are being pursued to 
measure the dynamics of M on small length scales. 
Among them, X-ray magnetic circular dichroism 
(XMCD) transmission microscopy^, XMCD photoelec- 
tron microscopy (PEEM) 3 -, microfocus Brillouin light 
scattering (BLS)^, time-resolved scanning Kerr mi- 
croscopy (TRSKM)£&£, spin-torque driven ferromag- 
netic resonance (ST-FMR)^! 1 ^ and magnetic resonance 
force microscopy (MRFM) applied to ferromagnetic 
resonanc o n i 12 ' 13 are the most accomplished yet. This 
article focuses on the last one, called herein mechanical- 
FMR, since a mechanical setup is used for the detection 
of FMR, as shown in FigQ] Here, we shall review the 
ability of mcchanical-FMR to extract quantitative spec- 
troscopic information in individual samples. 

In FMR, the magnitude of the magnetization vector, 
\M(t,r)\ = AI S , is a constant of the motion and equals 
the saturation magnetization at the lattice temperature. 
Thus the dynamics of M is governed by the position- and 
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FIG. 1: (Color online) Schematic representation of the 
mechanical-FMR setup: the resonance spectra of a small size 
sample are detected by a magnetic force microscope. The 
probe consists here of a magnetic sphere glued at the extrem- 
ity of the cantilever. The excitation comes from a microwave 
antenna placed underneath the sample. 



time-dependent transverse fluctuations of the magnetiza- 
tion direction, the so-called spin- waves (SW)^ m(t, r) = 
M{t,r)-M c (r)C,(r), where M c = M-£ is the longitudi- 
nal component of the magnetization, defined as the pro- 
jection of the instantaneous vector M(t) along the local 
precession (or equilibrium) axis of unit vector C( r )- Be- 
cause of the exchange interaction, these transverse fluc- 
tuations are a collective precession, and become more in- 
sightful once decomposed along the normal modes basis 
of the sample. Therefore the experimental identification 
of this SW basis is important. Its nature depends primar- 
ily on the symmetry of the equilibrium configuration of 
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M inside the sample^. In the case of a uniformly mag- 
netized sample, the SW confinement is mostly governed 
by the shape of the sample. In thin films, the quantiza- 
tion of the lowest energy (or longest wavelength) modes 
develops along the finite thickness direction, whereas in 
nanostructures, it arises both from the thickness and lat- 
eral confinements. 

Although, there arc fundamental interests in study- 
ing the magnetization dynamics in the time domain (or 
free precession regime) , the best spectral resolution is ob- 
tained by continuous wave spectrometers, which study 
the harmonic response (or forced oscillation regime). 
Conventional spectrometers excite and detect the SW 
cigen-modes of a ferromagnetic sample with the same mi- 
crowave antenna. In a conventional-FMR experiment^, 
the sample is placed in a region where a large homo- 
geneous static magnetic field aligns the spins in a well- 
defined direction. On the emission side, a small mi- 
crowave field h, applied perpendicularly to the effective 
magnetic field JT e fiF (|| C)> wu l cant the magnetization 
away from its equilibrium axis if the resonance condition 
is met, i.e. if the photon energy (fuj s ) corresponds to the 
energy to excite one SW or magnon (fvyH^s). The sym- 
metry of the excitation antenna sets the selection rules 
and the symmetry of the excited SW. On the reception 
side, the measured quantity in a conventional spectrom- 
eter is P a bs, the microwave energy absorbed inside the 
sample per unit of time. Technically, P a b s is obtained by 
monitoring the transmitted or reflected power to a mi- 
crowave diode, whose dc voltage is proportional to the 
incident microwave power. In metallic samples, the ab- 
sorbed power per unit of surface can be expressed as a 
function of Z s , the complex surface impedance of the 
sampl e 17 ' 18 , which itself depends on the complex elec- 
trical conductivity and magnetic susceptibility For thin 
magnetic layers, it mainly reduces to the dissipative part 
x" of the microwave transverse magnetic susceptibility, 
as in insulators: P a b s = J v dtM.h dV = uj s x"h 2 , where 
the shorthand notation dtM represents the time deriva- 
tive. 

The dependence of P a b s on cither the applied magnetic 
field, P a ppi, or the frequency of the excitation source, 
u) s /2n, reveals resonance peaks (bell-shaped curves). 
Their positions, amplitudes and line widths reflect the 
complete spectral information about the spin system, i) 
The position is a precise measurement of the effective 
magnetic field defined as the conjugate variable of the 
magnetization, H c g = 8mJ~, where T is the free en- 
ergy of the spin system, ii) The amplitude is related 
to the coupling (overlap integral) with both the exci- 
tation and the detection schemes. It thus gives a hint 
about the spatio-temporal profile of the mode. For in- 
stance, with a cavity (cavity-FMR) or a stripline antenna 
(striplinc-FMR), the microwave field h is uniform over 
the sample volume. As a result, it preferentially cou- 
ples to the longest wavelength SW modes of the sample, 
iii) The frequency line width gives the decay rate of the 
excited SW (coherent with h) to the other degrees of free- 



dom: the degenerate and thermal SW (incoherent with 
h) or the lattice (electrons and phonons). In Bloember- 
gen and Wang's notation^, the line width is proportional 
to Experimentally, it is often easier to monitor 

the dependence of P a b s when P a ppi is swept at fixed ui s . 
In this case, it is crucial to repeat the measurement at 
several frequencies and to renormalize the field line width 
by the effective gyromagnetic ratio 7 e fi = dui s / dH &vv \. 

FMR spectroscopy is usually restricted to probe tiny 
deviations of the magnetization from its equilibrium. 
First, reaching large angles of precession requires a lot 
of microwave energy, especially in metallic samples and 
second, the interpretation of the high power regime re- 
quires a specific analysis since the equation of motion is 
non-linear—. The typical angles of precession found in 
FMR are less than § = 1° . Although FMR is a sensitive 
technique, still the sensitivity of most conventional spec- 
trometers is usually not sufficient to detect magnetization 
dynamics in individual sub- micron size samples. The cru- 
cial parameter here is the ratio between the volume of the 
sample and the volume of the detector: the so-called fill- 
ing factor. With the recent development of MRFM, the 
sensitivity of magnetic resonance has been tremendously 
enhanced. Nanometer scale sizes come now within reach 
as Dan Rugar and colleagues at IBM Almaden recently 
showed with the detection of a single electron spirit and 
with a record 90 nm resolution on nuclear paramagncts, 
where only a few thousands of nuclear spins contribute 
to the signal^. The key to such performance comes from 
the size of the magnetic probe, which has a very strong 
coupling with the resonant volume, thus ensuring an al- 
most optimized filling factor. 

In this paper, we show how a mcchanical-FMR setup 
can be used to perform FMR spectroscopy of sub-micron 
size samples. We present a thorough study of the SW 
eigen-modes in individual permalloy (Py) disks, whose 
FMR spectra are compared to a 2D approximate analyt- 
ical model and to full 3D micromagnetic simulations of 
the dynamical susceptibility The structure of the rest 
of the paper is the following. Section [II] presents the 
principles of the mcchanical-FMR and gives an extensive 
analysis of the spectral deformations induced by the mag- 
netic probe on the measured FMR spectra. In scction HTTl 
the experimental results are presented. In section HVl we 
derive an approximate 2D analytical model and present 
the results of a 3D micromagnetic simulation. They are 
both used in section [V] to analyze the data. Section PVTl 
contains a summary of the results obtained. 



II. MECHANIC AL-FMR 
A. The experimental setup 

The concept of mechanical detection of the magnetic 
resonance was first applied to FMR by Phil Wigcn 
and Chris Hammcl in 1996ii. Figfl] illustrates the 
mechanical-FMR setup, while Figf2] presents its cxper- 
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FIG. 2: (Color online) a) Photography of the MRFM bottom 
base showing (1) the microwave antenna, (2) the mica sub- 
strate where the Py disks are, (3) the piezo-electric scanner, 
and (4) the magnetic coil used for a calibrated modulation of 
the applied field, b) Section view of the above elements. The 
top part with the cantilever and the optical detection is not 
shown. 

imental realization for the present study. 

1. Excitation part 

The sample magnetization is excited by a microstrip 
antenna placed underneath (see Figs Q] and The an- 
tenna consists of a 1 (im thick gold micro-strip patterned 
by optical lithography on top of a sapphire substrate 
whose bottom surface is a ground plane (i.e. fully cov- 
ered by Au) . In order to obtain a wide-band antenna, the 
top Au electrode is shorted to the ground plane at the 
extremity of the substrate, hereby creating an anti-node 
of microwave magnetic field at the extremity. Placing the 
sub-micron Py disks at this location ensures a wide-band 
excitation scheme, at least until the magnetic field node, 
a quarter of the wavelength away, moves underneath the 
samples, which occurs at very high frequency (above 20 
GHz). The efficiency of this wide- band setup to induce 
large variation of the absorbed power, depends i) on the 
amplitude h of the excitation, ii) on the volume of the 
sample, and iii) on the susceptibility x" of the sample at 
resonance, inversely proportional to the line width. It is 
important to keep the amplitude of h below saturation 
or other critical thresholds^. Indeed, significant spectral 
distortions arc induced due to non- linear effects. These 
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FIG. 3: (Color online) Visualization of the dipolar coupling 
interaction between a spherical magnetic probe (diameter ^) 
and a cylindrical Py disk (diameter D) separated by a distance 
s. The magnetic configurations in both objects are assumed 
to be saturated along z. The color lines represent the iso--B z 
field lines. The bright and dark shades separate respectively 
the positive and negative coupling regions. 



distortions being themselves a subject of research^, all 
the data shown hereafter arc taken in the linear regime, 
where the peak amplitude remains proportional to the 
excitation power i.e. precession angles limited to 1°. 



2. Detection part 

The detection scheme is directly inspired from mag- 
netic force microscopy (MFM). 

A soft cantilever with a magnetic tip is placed in the 
stray field of the sample to be studied. The dipolar cou- 
pling between the tip and the sample (see Fig(3j) creates 
a small flexion of the cantilever beam, which is detected 
optically by deflection of a laser beam on a 4-quadrants 
photodiode. The pitch angle is produced by the vertical 
component of the force F z (see FigQ} but also by the 
torque N y , the component perpendicular to the beam 
axis and to the vertical direction. In all the experiments 
below, we will measure with a spherical probe the mag- 
netization dynamics of disk-shaped Py samples perpen- 
dicularly magnetized by an homogeneous applied field. 
H&ppi is larger than the saturation field of Py, such that 
the magnetic state of the sample is almost uniform. In 
this situation, the coupling geometry preserves the ax- 
ial symmetry, hence the mechanical coupling comes from 
the force F z only (N y = 0). 

There is a principal difference between the detected 
signal obtained by mechanical-FMR and by conventional- 
FMR. Mechanical-FMR uses detectors which are only 
sensitive to M^, the static (or longitudinal) component 
of the magnetization, i.e. to the susceptibility averaged 
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on time scales which are much longer than any relax- 
ation time in the spin system. The microwave oscilla- 
tions of the transverse component do not couple to any 
mechanical mode of the cantilever, which typically oscil- 
lates in the audio range. When an FMR resonance is 
excited by the antenna, it is the deviation of Mq from 
M s , AM C = M s - M ( , that diminishes the force F z 
on the cantilever (see FigfT]). In a corpuscular picture, 
each magnon reduces M^ by 7/1. Therefore the variation 
AM^ is proportional to the total number of magnons 
being excited by the microwave field inside the sample, 
independently of their coherence. Another way to ex- 
press this is through energy arguments. While conven- 
tional FMR measures the absorbed power (see above), 
mechanical-FMR measures the stored energy^. The ra- 
tio between these two quantities is actually 1/Ti, the re- 
laxation rate of the out-of-equilibrium full magnons' pop- 
ulation towards the lattice. In other words, mechanical- 
FMR provides an intrinsic information about the dynam- 
ics of the SW system, where the area under each reso- 
nance peak of the AM^ spectrum is proportional to T x . 
But the AMq measurement is seldom obtained, because 
this quantity is a second order effect in the precession 
angle AMq/M s = 1 — cosiS is thus much smaller than 
the transverse susceptibility, which is of first order in 



3. Imaging 

Quite naturally, the scanning probe can be used here 
for imaging purpose. But the situation is different in 
paramagnetic^ and ferromagnetic spin systems. In a 
paramagnet, the spins are decoupled and the excitation 
is localized in a "resonant slice"—, whose thickness is 
inversely proportional to the field gradient produced by 
the magnetic tip. In contrast, neighboring spins in a fer- 
romagnet are coupled through the exchange interaction, 
which works against localization effects. The stray field 
of the magnetic tip can, however, alter the sample mag- 
netic configuration just underneath the tip, leading to 
tip-induced FMR resonances. The group of Chris Ham- 
mel has recently showed, that it is possible to produce 
these new localized FMR modes by approaching the mag- 
netic tip close to the sample surface^. Exploitation of 
the spectral features of these new modes is still a chal- 
lenge at the moment. 

Detection wise, one benefits from the same advantage 
of an MFM. The spatial resolution is related to the size of 
the magnetic probe and the separation with the sample. 
In the present study, we are mostly interested by an op- 
timization of the sensitivity of the mechanical detection, 
allowing us to measure smaller sample sizes. As will be 
shown below, optimization of the sensitivity requires to 
chose a size of the magnetic probe of the order of the size 
of the sample (optimum filling factor—). Furthermore, 
placing the probe far away from the sample surface and 
working in the weak coupling regime diminishes the spec- 
tral alteration produced by the tip. Such conditions are 



quite obviously incompatible with good imaging condi- 
tions and an increase of the spatial resolution (smaller 
probe) must then come at the detriment of the sensitiv- 
ity. 



B. Detection sensitivity 

One of the main advantage of the mechanical-FMR 
is its exquisite sensitivity. This is mainly due to the 
progress in nanofabrication technologies, which can pro- 
duce micron-size mechanical structures with outstanding 
performance figures, i.e. achieving among the best com- 
promise between small size and large quality factor. 



1. Modulation technique 

Exciting the sample at a fixed frequency (u) s /2ir), spec- 
troscopy is achieved by recording the cantilever motion 
as a function of the perpendicular dc applied field, -ff a ppi, 
produced by an electromagnet. As mentioned before, the 
coupling between the mechanical oscillator and the mi- 
crowave magnetization dynamics is purely static. How- 
ever, it is possible and useful to modulate the microwave 
power at the mechanical resonance frequency of the first 
flcxural mode of the cantilever f c = 3 kHz in our case). 
As a result, the amplitude of vibration will be multi- 
plied by Q, the quality factor of the mechanical resonator. 
Note that the mechanical noise is also amplified. How- 
ever, a sensitivity gain is obtained if this intrinsic me- 
chanical noise exceeds the preamplifier noise or the noise 
of the microwave source. This modulation technique is 
referred to source or amplitude modulation. The ampli- 
tude of the microwave field follows the time dependence 

ft(t) = fte <w '*|i + ico 8 (27r/ c t)} ) (1) 

where the depth of the modulation is 100%. Note that 
this modulation technique does not affect the line shape 
in the linear regime, because the period of modulation 
l// c , is very large compared to the relaxation times T\ 
and Ti of the ferromagnetic system studied. 

We mention that the modulation of the microwave field 
at f c also induces a direct vibration of the cantilever, even 
outside any resonance phenomena of the probe or of the 
sample. We attribute this to a modulation of the temper- 
ature of the cantilever. The latter is a direct consequence 
of the modulation of the microwave heating and eddy cur- 
rents, mainly induced by e, the electric component of the 
electromagnetic radiation. This effect can distort the res- 
onance peaks due to the modulation of the applied field 
on the sample by the vibrating magnetic probe. For in- 
stance, a sphere with a magnetic moment of 6 10 -9 emu 
vibrating by 10 nm pp induces a field modulation of about 
6 Ocpp on a sample placed at 2.7 /im from its center (sep- 
aration s = 1.0 /im, see below). This field modulation 
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has no influence on the FMR signal if it is small com- 
pared to the line width, but will otherwise significantly 
broaden the resonance. To eliminate it, a forced oscilla- 
tion out of phase with the spurious contribution is pro- 
duced by a piezo-electric bimorph slab placed nearby the 
cantilever. The resulting total vibration corresponds to 
a few nanometers at most. All the spectra shown below 
are obtained in these conditions, where the amplitude of 
vibration of the cantilever is compensated so that it has 
no influence on the measured FMR spectra. 

2. Minimal detectable force 

In all MRFM setups, the detection noise is only limited 
by the Brownian motion of the cantilever, which behaves 
as an harmonic oscillator with a single degree of freedom. 
The theorem of equipartition of energy stipulates that, 
in the absence of external force, the vibration amplitude 
of the cantilever is such that the average kinetic energy 
associated to the thermal excitation is equal to ksT/2, 
where fcs is the Boltzmann constant and T is the tem- 
perature of the bath. Thus the minimal detectable force 
follows the relation 

/ 2k B TkB 

rmin — \ , „ , (2) 

where k is the spring constant of the cantilever, Q its 
quality factor, and B the detection bandwidth. Fem- 
toNewton sensitive cantilevers are now readily available 
commercially. A BioLcvcr B from Olympus with k = 
5 10~ 3 N/m has been used for this work. Because the 
cantilever is very sensitive to thermal fluctuations, it is 
important to stabilize both the intensity of the laser (to 
less than 15 ppm) that is shined on top of the cantilever 
for the position sensing but also the over-all temperature 
of the microscope (to less than 200 ppm) , that is mounted 
on Peltier elements. To increase the sensitivity, the ex- 
periment is operated in a secondary vacuum of about 
10~ 6 torr. Operating in vacuum is important to keep the 
large value of Q = 4500. As seen in Eq.©, the other 
parameters that control the minimum detectable force 
are the temperature and the bandwidth. Working at low 
temperature and averaging the signal over large period 
of time allows the detection of attoNewton forces^ 9 -. All 
the data shown hereafter are obtained at T = 280 K and 
with a lock-in time constant of one second. 



3. Magnetic spherical probe 

Obviously the sensitivity directly depends on the 
strength of the dipolar coupling between the probe and 
the sample. As mentioned above, the magnetic probe 
must be carefully chosen, because there is an optimal 
size for a given sample (concept of filling factor). In our 
case, the probe itself is a magnetic sphere glued at the 
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FIG. 4: (Color online) Variation of the dipolar magnetic force 
between a sphere and a disk in contact (s = 0) as a function 
of the ratio between their diameters, respectively ^ and D. 
The maximum force occurs when = 2D. 

apex of the cantilever. For a spherical probe, the static 
force applied on the cantilever is simply proportional to 
g zz , the z component of the field gradient created by the 
sample at the center of the sphere: 

F z = m sph gzz, (3) 

where m sp h is the magnetic moment of the sphere. For 
a disk-shaped sample, the largest force is obtained when 
the radius of the sphere is equal to the diameter of the 
disk^S, as shown in Fig[4j The optimum corresponds to 
the particular case where the sphere captures all the pos- 
itive field lines emanating from the disk (see Fig [3]). 

A Scanning Electron Microscopy (SEM) image of the 
sphere is shown in Figj5] Its diameter is ^ = 3.5 /im, 
which is slightly larger than the optimal size for our 
D\ = 1.0 and D2 = 0.5 /itm sample diameters. This 
sphere is an amorphous alloy whose main constituents 
arc Co (80 wt%), Fe (10 wt%) and Si (9 wt%), as de- 
duced from chemical analysis. Its characteristics are 
measured after the gluing process of the sphere at the 
tip of the cantilever. The magnetization curve is ob- 
tained by placing the mounted cantilever above an Fe 
cylinder (diameter 2 mm, height 8 mm), which creates a 
well-characterized field gradient of g zz = 0.5 G/fim. By 
monitoring the deformation of the cantilever vs. the ap- 
plied field, we can infer a 1.2 kG saturation field for the 
sphere. We also obtain the value of its magnetic moment 
"^sph = (6 ± 0.5) 10 -9 emu, which is deduced from the 
known values of the field gradient and the cantilever's 
spring constant. 

We have already mentioned above the interest to pre- 
serve the axial symmetry, which is achieved if the center 
of the sphere is positioned above the center of the disk. 
Two other comments should be added about the choice 
of the spherical shape for the magnetic probe. First, 
in contrast to cylinders, spheres do not have any shape 
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FIG. 5: (Color online) Scanning Electron Microscopy image 
of the CoFeSi magnetic sphere glued at the tip of an Olympus 
Biolever. 

anisotropy. This property, combined with the absence of 
magneto-crystalline anisotropy of the amorphous alloy, 
considerably decreases the influence of the applied mag- 
netic field on the mechanical resonance frequency of the 
cantilever. This is important because the source modu- 
lation technique relies on a phase locked loop (PLL) to 
keep the source modulation exactly at the resonance fre- 
quency of the cantilever. Second, the magnetic probe on 
the cantilever is itself placed in the field gradient of the 
sample. As a result, mechanical vibrations of the can- 
tilever are also produced by the excitation of the FMR 
modes of the probe itself. This can be seen as a recipro- 
cal effect of the probe-sample coupling. Using different 
shapes for the probe (sphere) and for the sample (disk), 
the two FMR spectra are completely separated, mainly 
due to the difference in the shape anisotropics of the two 
different geometries. In the spectral range shown in this 
study (w s /27r = 4 to 10 GHz), the FMR resonance of 
the sphere occurs below 3 kOe 3 -^, i.e. at much lower field 
than for the disks studied hereafter. 



C. Extrinsic effects 

1. Spectral deformations induced by the magnetic probe 

The optimization of the sensitivity discussed above 
has, however, a drawback. Increasing the coupling be- 
tween the probe and the sample inevitably produces a 
distortion of the spectral features. It is thus important 
to understand this effect and to keep it inside the per- 
turbation regime so that one could still deduce the in- 
trinsic behavior of the sample, i.e. without the presence 
of the probe. In mechanical-FMR, the spectral distor- 
tions are due to the stray field of the probe. It produces 
an additional field inhomogeneity in the internal effective 
magnetic field, which affects the detailed SW dispersion 
relation inside the sample and thus the resonance condi- 
tion. The strength of the effect depends on m sp h and on 
the separation s between the probe and the sample (see 
Figs U and [3]). We discuss in the following how to en- 



sure that the field inhomogeneity produced by the probe 
is small compared to the internal dipolar field variations 
along the radial direction of the sample (« 2irM s for a 
non-ellipsoidal shape). 

Thorough experimental studies have been performed in 
the past on mechanical-FMR of a YaFesO]^ single crystal 
disk for different separations s between the magnetic tip 
and sample surface^. It was found that, when the bias 
field inhomogeneity from the probe is smaller than a few 
pcrcents of the internal field variation inside the sample, 
its main effect is to shift the entire spectrum to higher fre- 
quency as s decreases. This shift is homogeneous within 
less than 10% for all the peaks corresponding to the dif- 
ferent SW modes. In these conditions, their spatial pro- 
files are almost not affected by the presence of the probe, 
and the relative amplitudes between the peaks are kept 
to their intrinsic values. On the other hand, if the bias 
field from the tip represents a large perturbation to the 
internal field inside the isolated sample, the spatial pro- 
files of the SW modes are affecte d 31 ' 32 , and the observed 
FMR spectrum can not be recognized as intrinsic to the 
sample. This situation can also lead to a different cou- 
pling between the probe and the sample 3 -^. However, it 
is to note that since the mechanical probe is independent 
from the excitation, an increase of the coupling between 
the sample and the probe do not produce any additional 
contribution to the FMR line widtb-^. This is in contrast 
to the radiation damping found when the coupling to the 
microwave resonator increases^. 



2. Quantitative analysis 

A quantitative estimation of -ff ff, the shift induced by 
the tip on each SW mode, can be obtained analytically 
within a 2D model. By 2D, we mean here a model where 
the precessional profile of each mode in the disk only de- 
pends on the two in-plane coordinates, i.e. on (r,cf>) in 
a cylindrical frame (it is uniform along the thickness di- 
rection z). In that case, the resonance condition depends 
on the phase delay of the SW accumulated along the di- 
ameter. Constructive interferences occur when the phase 
delay over a cycle is equal to n x 2tt, n £ N*. This condi- 
tion is equivalent to the WKB approximation applied to 
the dispersion relation of magnetostatic forward volume 
waves (MSFVW), established by Damon and Eshbach in 
1961—. In the following, we will use the dipole-exchangc 
dispersion relation developed by Kalinikos and Slaving 
for magnetized thin films and later applied by Kakazei et 
al. for the case of disks magnetized in the exact perpen- 
dicular geometry 3 -^. 

We start with the normal modes basis of a disk mag- 
netized in the perpendicular direction, which are the 
J^m( r ) = J^{ki, m r) cos(£(/)), where Jgs are the Bessel func- 
tions of the first kind and kg^ m is the modulus of the 
in-plane SW wave vector determined by the boundary 
conditions. We shall assume that the modes satisfy the 
dipolar pinning condition at the circumference of the disk 
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FIG. 6: (Color online) Color code representation of the eigen- 
modes basis Sm(r), where I and m indicate respectively the 
number of nodes in the circumferential and radial direction. 
Only the modes 1 = couple to a uniform microwave excita- 
tion. 



samples, i.e. J'^(R) = 0. Thus fc^ m = Ke yTn /R, where 
K£. m is the (m + l) th root of Je(x) and R is the radius of 
the disk. In this notation, £ and m are respectively the 
azimuthal and radial mode indices (i.e. the number of 
nodes in the circumferential and radial directions). Fig[5] 
is a color-coded representation of the transverse suscep- 
tibility x" corresponding to the first modes. 

The obtention of an analytical formula for the reso- 
nance condition shall be detailed in section IIV Al for the 
general case where the applied magnetic field makes a 
small angle Oh with the normal of the disk (see Eqs. (f29|) 
and (|5CT|) ). For the estimation of the spectral defor- 
mations induced by the probe, we are only interested 
by the Oh = 0° limit of the linearization of the equa- 
tion of motion. The resonance condition for the normal 
mode (£, m) in the perpendicular direction is given by the 
expression^: 



w s 



, W int + W A /A 2 fc2 



+ uj m A.% m + - GtJ^M ) ■ (4) 

Here, cum = 4:irjM s and Wi n t = jH mt is the mode depen- 
dent effective internal magnetic field 



H A , 



(5) 



where H cxt is the total external magnetic field, N zz is 
the longitudinal matrix element of the effective demag- 
netizing tensor (see Eq.©) and Ha is the perpendicular 
uniaxial anisotropy field (of spin-orbit coupling origin). 



A = \J2A/ '(AttM^) is the exchange length, which de- 
pends on the exchange stiffness constant A, expressed in 
erg/cm (= 10~ 6 in Py). Finally, the parameter G^ m is 
derived in section [TV Al Eq. (f!?rj|) . 

In our notation, N is the demagnetizing tensor of the 
disk. Appendix [Al gives the analytical expression for the 
different matrix elements of this tensor in the cylindrical 
coordinates. In Eq.([5|), the demagnetizing field along z 
depends on a matrix element of the effective demagne- 
tization tensor iV^. m . The latter is the demagnetizing 
tensor weighted by the spatial dependence of the normal 
mode profile jf n {r), 



N, 



1 



W,m J r <R 

where Ce m is a renormalization constant 



iV(r)j£(r-)Vr, (6) 



'-'rr 



) 2 d 2 1 



(7) 



r<R 



To calculate the influence of the sphere on the reso- 
nance field of the mode (£, m), we need to expand the 
external field as the sum of two contributions: the homo- 
geneous magnetic field H app \z produced by the electro- 
magnet and the inhomogeneous stray field of the sphere 
Hsph(r): 



H eKt (r) = .ffappi + H sph (r). 



(8) 



Here we are only interested by its z-component, 
H sp h(r) = H sp h(r) ■ z. Then, the influence of the probe 
simply yields a modification of the internal magnetic field 
Hint of Eq.© by 



Hint = 



appl 



H oS - N zz 4wM s + H A , 



(9) 



where the shift H s induced by the probe is averaged 
along the radial direction by the mode profile: 



H 



1 



off 



C, 



r<R 



H sph (r)ji(rfd 2 r. 



(10) 



This expression allows us to give an estimation of the 
maximum coupling allowed to keep the spectral defor- 
mations in the perturbation regime. The additional con- 
tribution H Q g must be small compared to the internal 
field variation inside the sample volume (~ N zz 4:irM s ) in 
order to keep the normal mode basis unchanged. In the 
opposite case, the SW mode profiles and their resonance 
fields must be directly calculated from the MSFVW dis- 
persion relation and from the magnetostatic potential. 
Clear experimental signatures of this strong influence of 
the probe on the SW mode can be found^L. The experi- 
ments presented below are not in this regime. 

Quantitative results of the calculation are displayed in 
Fig|7]for the two different disks diameters, D\ = 1.0 and 
D 2 = 0.5 fxm. We have calculated the influence of the 
probe on the circumfcrentially symmetric (£ = 0) normal 
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FIG. 7: (Color online) Calculated spectral deformation in- 
duced by the magnetic tip placed at distances s = 3.1 and 
1.0 /im above the disk samples D± — 1.0 and D2 — 0.5 fim, 
respectively. The graphs show the downward shift in field of 
each resonance mode (indexed by m) compared to the unper- 
turbed FMR spectrum (dashed line) . 



mode J^{r). The value of M s and 7 used in the calcu- 
lation are those derived below from cavity-FMR studies 
of the extended thin film. We also use the physical char- 
acteristics of the probe derived from the SEM images in 
FigJS]and magnetometry. The values of s (the separation 
between the sphere and the disk) used in Fig[7] are close 
to the experimental conditions. The magnetic probe is 
brought closer to the disk surface as the sample diameter 
diminishes in an attempt to compensate for the reduction 
in the coupling strength. The induced force is expected 
to decrease by approximately an order of magnitude be- 
tween the two disks diameters: a factor of 4 comes from 
the decrease of the resonating volume and a factor of 2 
comes from the decrease of the filling factor (see dots on 
FigH]). On the other hand, approaching the probe from 
s = 3.1 to 1.0 fj,m corresponds to an increase of the gra- 
dient by an order of magnitude (from g zz = 6.1 x 10 5 
to 5.8 x 10 6 G/cm, respectively). For comparison, the 
dashed line is the s = 100 fim case, which represents 
the intrinsic, unperturbed case. Two important observa- 
tions can be deduced from these plots. First, the shift is 
almost independent of the mode number m. The spec- 
tral deformation can thus simply be modeled by an offset 
on the applied field. Second, at constant coupling, the 
amplitude of this offset increases when the disk diame- 
ter decreases, as the probe has to be brought closer to 
the smallest sample. In the case of our magnetic sphere, 
the displacement of the resonance field is about —100 Oc 
downward in field for the large disk and of —520 Oe for 
the small one. 

In conclusion, under experimental conditions of an al- 
most optimized coupling between the sample and the 
probe, the influence of the latter on the intrinsic FMR 
spectrum of the former is only an over-all downward shift 



in field, which can be quantitatively estimated. The fact 
that our sphere has a small moment and a large diame- 
ter keep the stray homogeneous over the sample volume 
and reduces deformation in the relative position of the 
resonance modes. Assuming that the sphere is brought 
into contact with the surface of the disk, the variation of 
the perpendicular stray field along the radial direction of 
the disk is about 400 Oc, which is small compared to the 
N zz AirM s ss 10 4 Oe variation of the internal field inside 
the disk. 



III. SPECTROSCOPY OF SUB-MICRON SIZE 
DISKS 



A. Film layered structure 



Films corresponding to the composite system (permal- 
loy = Ni8oFe2o, abbreviated by Py | alumina) have been 
deposited 3 ^ by RF sputtering, at room temperature, on 
two different substrates: single crystalline Si and mica. 
A sweeping mode for sputtering has been used, where 
for each material the substrate-holder is swept back and 
forth over the activated target. Such a mode of deposi- 
tion was chosen to achieve a better homogeneity of the 
magnetic layer and a good control of its thickness. Sev- 
eral multilayers (Py | A^Oa)^ have been produced, us- 
ing this sweeping mode, being characterized by an in- 
dividual layer in the ultrathin range. Low-angle X-ray 
diffraction measurements have been performed on these 
multilayers (Si substrate) and the exploitation of these 
results has led to the determination of the deposition 
rates for the two materials. The specific magnetic film 
studied in the present report has the following layered 
structure: (AI2O3 base | Py | AI2O3 top). It consists of a 
single Py layer (43.3 nm thick) sandwiched by two AI2O3 
layers of identical thickness (16 nm). The top alumina 
layer protects the Py layer from oxydation. The smooth 
amorphous AI2O3 base layer, characterized by a small 
surface roughness (a few A) helps the growth of (111) 
textured Py polycrystallinc layers in the low thickness 
range for the Py layers. The typical spread of surface 
orientation of the crystallites in the (111) textured Py 
layer grown by sputtering on the amorphous alumina is 
A9h = 0.5°. Due to the symmetry of the multilayer 
structure, the pinnings of the magnetization at the top 
and bottom (Py | AI2O3) interfaces are expected to be 
identical. Most importantly it is to be noted that, in this 
layered structure, the Py layer is sandwiched by an in- 
sulator. Consequently, by contrast to the cas o 37 ' 38 of an 
adjacent metallic layer (Py | N), where N designates a 
normal metal, no modification of the Py intrinsic damp- 
ing (a spin-diffusion phenomena) is expected to arise from 
the (Py I AI2O3) interfaces. 



B. Cavity-FMR studies of the extended thin films 

Cavity-FMR experiments have been performed on the 
extended thin films deposited on Si and mica, corre- 
sponding to the layered structure described above (16 
A1 2 3 I 43.3 Py I 16 AI2O3), with individual thickness in 
nm. The basic FMR experiment which has been carried 
out (reflexion at X-band, 9.6 GHz, and room temper- 
ature) consists in studying the resonance spectrum as a 
function of the orientation Oh = (z, H app \) of the dc field 
-ffappi applied in a plane perpendicular to the film. The 
resonance condition^ for the uniform mode, as a function 
of 9h, depends only on two parameters: the gyromag- 
netic ratio 7 and the total perpendicular anisotropy field, 
H u = 4ttM s —Ha- The two fitting parameters (7, H u ) are 
deduced from the observed angular variation of the reso- 
nance field. Very close values of these two parameters are 
found for the the films grown on Si and on mica. For the 
film grown on mica, used to pattern the sub-micron disks, 
the gyromagnetic ratio is 7 = 1.849 10 7 rad.scc _1 .Oe _1 , 
corresponding to a Lande factor g = 2.103 ± 0.004. The 
perpendicular anisotropy field H u = 9775 ±50 Oe reflects 
entirely the demagnetizing field 4nM s and corresponds to 
the expected value of the magnetization for a NiFe alloy 
of atomic composition NigoFc2o. In fact, the spin-orbit 
anisotropy field Ha, that is ascribed^ to a stress- induced 
anisotropy observed in the ultra-thin NiFe layers, is here, 
for this thick Py layer, nearly zero (at least < 100 Oc) 
due to the fact that for this alloy composition the mag- 
netostriction reduces nearly to zero. We mention also 
that this thick Py layer is characterized by a small uni- 
axial in-plane anisotropy, which is identified to a field- 
induced anisotropy built during the film deposition. Its 
small value, 7.5 ± 0.5 Oe, is deduced from the two paral- 
lel geometry FMR measurements with the magnetic field 
applied respectively along the easy and the hard axis. 



C. Mechanical-FMR studies of the patterned 
structures 

The disks prepared for the mcchanical-FMR experi- 
ment were patterned out from the (AI2O3 | Py | AI2O3) 
film grown on mica. Because this substrate cleaves eas- 
ily, it was possible to reduce its thickness down to about 
15 //m before gluing it on the broadband stripline. As 
a result from the short distance between the sample and 
the excitation circuit, the ac field h at the sample location 
can achieve up to 10 Oe for the available power from the 
synthesizer. However, in the results presented below, h 
was restricted to 1 Oe in order to avoid non-linear effects. 
The mcchanical-FMR data are all collected at T = 280 
K. The microscope was aligned so that the normal of the 
mica substrate was lying within 5° with the dc field. 

The disks studied below have been patterned out of 
the same thin film deposited on mica which was studied 
by cavity-FMR. An aluminium mask defined by e-beam 
lithography is used to protect the Py disks during the 
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FIG. 8: (Color online) SEM images of the two Py disks sam- 
ples. The images are being viewed at an angle, thus distor- 
sions exist in the vertical direction. 



subsequent ion etching of the thin magnetic film. Several 
disks of nominal diameters 0.5 and 1.0 /xm separated from 
each other by 50 [im. were defined using this lithographic 
process. Fig|S] shows SEM images of two such patterned 
disks. They do not correspond to the measured disks, 
whose FMR spectra are presented in Figsl9l andfTTJl but 
were processed using the same recipe. A thin Au layer 
had to be deposited before the SEM imaging to avoid 
charging too much the mica substrate. The shape and the 
dimensions of the samples are as expected (the slightly 
elliptical shape seen on the images is due to perspective) . 
The lateral rugosity, which is about 30 nm, could be due 
to uncertainties in the lithographic process on mica, or 
to an inappropriate dose during the insulation. It is not 
expected to influence too much the eigen-modes profiles 
of a perfect disk (presented in Figj6j) and their resonance 
fields as their amplitude vanish at the periphery of the 
disk. These radial fluctuations do not either contribute 
to the FMR line width when the measurement is made 
with the field parallel to the normal of the disk. In that 
case, only inhomogencitics along the thickness (i.e. the 
magnetization direction) have a direct impact on the line 
width broadening^. 
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FIG. 9: (Color online) Mechanical-FMR spectra of the Di = 

1.0 /im disk measured at 4.2, 5.6, 7.0 and 8.2 GHz and s ~ 

3.1 /im. The lines are the analytical predictions of the locus 
of the (l, m)-modes when 9h = 4.9° and H s = —70 Oe. 



FIG. 10: (Color online) Mechanical-FMR spectra of the D 2 = 
0.5 fim disk measured at 4.2, 5.6, 7.0 and 8.2 GHz and s ~ 
1.0 /im. The lines are the analytical predictions of the locus 
of the (I, m)-modes when 6h = 9.5° and H r = —760 Oe. 



1. Mechanical-FMR study of the D 1 = 1.0 \im disk 

The mechanical-FMR measurement on the D\ = 
1.0 /im disk was performed with a probe-sample sepa- 
ration s w 3.1 /im, kept constant at the four different 
frequencies studied and shown in Fig|9] To ensure that 
the separation is the same between the different spec- 
tra, a fine tuning of s is operated at the beginning of 
each scan (same H app \) so as to keep the frequency of 
the cantilever identical. A series of magnetostatic modes 
is observed at each microwave frequency, with the most 
intense at the highest field, i.e. at the lowest energy. It 
corresponds to the mode whose coupling with the exci- 
tation is maximum, because it is the longest wavelength 
mode. We note that the 5.6 GHz spectrum has been mea- 
sured for a larger integration time than the other spectra, 
which explains its better signal to noise ratio and allows 
to clearly observe the small amplitude low field modes. 
The detailed analysis of the spectral features (position, 
amplitude and line width) of the D\ = 1.0 fim disk is 
done in section [Vj based on the models developed in sec- 
tion [TV] 



2. Mechanical-FMR study of the D2 = 0.5 fim disk 

The mechanical-FMR study on the D2 = 0.5 /im 
disk was performed with a probe-sample separation s w 
1.0 /im, kept constant at the four different frequencies 
studied and shown in Fig[l0l As for the D\ = 1.0 /im 
disk, different magnetostatic modes can be observed, 
with the most intense at high field. The main difference 
with the larger disk is that this largest peak happens at 
much lower field 2.5 kOc less). It can not be explained 
solely by the larger stray field from the probe, which is 



closer from the sample. Actually, most of this shift is 
due to finite size effects. The latter are also responsible 
for the larger separation between modes on the smaller 
diski^. A detailed analysis of the spectral features (po- 
sition, amplitude and line width) of this disk is done in 
section [V] 

We note in passing that the signal to noise ratio ob- 
tained at T — 280 K with a time constant of one second 
on this tiny sample is very decent. One has to remem- 
ber, that the angle of precession is less than 1°. This 
translates in a spin sensitivity of about 1000 spins for 
our mechanical-FMR setup. We also recall that for this 
sample, the probe is not optimal, since ^/D 2 = 7 > 2. 

IV. MODELING OF THE SPECTRA 

We shall present below two different approaches to ana- 
lyze the experimental data. The first one is a 2D approx- 
imate analytical model, which assumes an homogeneous 
magnetization dynamics along the thickness. The second 
one is a full 3D micromagnetic simulation of the dynami- 
cal susceptibility using a discretized representation of the 
0.5 /im disk with regular cubic cells of 3.9 nm lateral size 
and assuming an homogeneous magnetization dynamics 
inside each cell. 

The motivation is to develop a comprehensive frame- 
work to anayze the spectroscopic features observed ex- 
perimentally by mechanical-FMR. The parameters intro- 
duced in the models are the magnetic characteristics of 
the material as extracted from cavity-FMR of the thin 
film and introducing finite size effects (with the assumed 
diameters and thickness of the disks) and a small mis- 
alignment Or 7^ 0° between the applied dc field and the 
normal of the disks to fit the mechanical-FMR data. To 
predict the effects of the misalignment, we shall use a 
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FIG. 11: (Color online) Graphic representation of the Carte- 
sian frame (£, £) (blue) rotated in the direction of the equi- 
librium magnetization, when the applied field makes an angle 
Oh with the normal of the disk (z-axis). 



perturbation approach, assuming that the profile of the 
modes is unchanged compared to the perfect perpendic- 
ular alignment. By comparing the analytical models and 
the simulations, we will evaluate the range of validity of 
this approximation. 



A. 2D analytical model 

In this section, we derive analytically the resonance 
condition that applies to the case where the external 
magnetic field makes a small angle Oh with the normal 
of the disk, which breaks the axial symmetry. The nota- 
tions used below are defined in Figllll 

The idea is to perform a linearization of the equation 
of motion for the magnetization around the equilibrium 
configuration using the SW dispersion relation in the 
dipole-exchange regime. We start with the gyromagnetic 
term of the Landau-Lifshitz equation of motion for the 
magnetization: 



d t M = 1 [H cS x M] . 
The effective magnetic field 

H g — Hext + -Hex + -f^dip 



(11) 



(12) 



is the sum of three fields: i) H cxt , the external mag- 
netic field, which includes the stray field produced by 
the probe; ii) the exchange field 



H cx = 4ttA 2 V 2 M 



(13) 



A being the exchange length; iii) Jfdip the internal dipo- 
lar field of the sample: 



H dip (t,r) 



-4tt / G(r 

lr'<R 



r') ■ M(t,r')d 2 r'. (14) 



Here R is the radius of the disk, and G(r) is the dipolar 
Green tensor, that has the form 



G(r) = / G k e 



, d 2 k 
(27T) 2 



(15) 



which introduces non-local interactions. The Fourier 
image Gk of the dipolar Green tensor for the lowest 
SW branch (with uniform magnetization across the film 
thickness) is given by^i 



G k = P k zz + (1 - P k 



kk 

) ~k 2 ~ : 



where 



-ki 



Pk = 



kt 



(16) 



(17) 



z is the unit vector orthogonal to the disk plane, t is the 
thickness of the ferromagnetic film, and fe the SW wave 
vector. 

For a thin magnetic disk (with aspect ratio t/R < 1) 
we can neglect non-uniformities of the equilibrium mag- 
netization distribution M(r). It will be assumed that 
M(r) = MsC, where £ is a constant unit vector along 
the effective magnetic field direction. At equilibrium, a 
uniform magnetization will create a non-uniform dipolar 
field — 4irM s N ■ where AT is the position-dependent 
static demagnetization tensor defined in appendix [X] 
The unit vector £ is determined by the condition: 



H cxt - 4ttM s N ■ C 
This can be rewritten as: 



(18) 



2nM s N xx sin 29 - 2ttM s N zz sin 26 = 0, 

(19) 

allowing to determine the equilibrium angle of the mag- 
netization with respect to the normal of the disk. N is 
the averaged demagnetization tensor over the volume of 
the disk: 



N 



1 



irR 2 



N(r)d 2 r. 



(20) 



r<R 



For small perturbations around the equilibrium state, 
the transverse part of the magnetization m(t, r) = 
M(t, r) — M S C obeys the linear equation: 



d t m 

where 
h dip (t,r) 
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4ttM s N ■ n X m 



+ 1 M S [(h dip + 47rA 2 V 2 m) x C] , (21) 



-4tt / G(r 

lr'<R 



r')-m(t,r')d 2 r'. (22) 



We recall that for axially symmetric samples with neg- 
ligible thickness, the normal modes are of the form 



m^ m {t,r) = Jm{r)ne, m (t) = M k e, m r) cos 



.(*), 
(23) 
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where H£ tm is a left circularly polarized unit vector ro- 
tating at u) s , t and m are, respectively, the azimuthal 
and radial mode indices, Jg(x) are Bessel functions, and 
the wave numbers ki yTn are determined by the boundary 
conditions. 

For a thin disk one can assume that the mode profiles 
remain the same as in the limit t — > 0. Substituting 
Eq.([23} into Ea.([2T]). multiplying by Je(k e , m r) cos(£0), 
and averaging over the disk area, one obtains the usual 
differential equation for the unit vector \xi_ m : 
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where Ne. m has been defined in Eq. 



H cxt - 4^M a N^ m ■ Cj X n t>rn 
47T7M S ( G^ m ■ fie tm + (Aki tm ) 2 fie, m ) x(§4) The tensor G^ m is equal to 



where 



G, 



1 



,m J r <R Jr'<R 

Cf m zz + (1 — Cf_ rn ){xx) 



G(r - r')Je(k e , m r)J e (k(, m r')e 



fr'd 2 r 



Go — 



2nR 2 



Ci 



m JO 



kJ t 



'i-i{kR)Jt(k^ m R) - k^ m Jt-i(ki>, m R)J t (kR) \ 



U2 _ 1,2 



Pkkdk. 



(25) 



(26) 



Pfc has been defined in Eq. (|17p and the expression for numerically. The resonance field is obtained by solving 
the normalization constant C^ m is given in Eq.([7]) and Eqs. ([24|) (see appendix [Bj) . It yields 
simplifies to 

Ct, m = ttR 2 (Jt+i(kt, m R)Jt-i(k e , m R)) ■ (27) 



Eqs. (f2"4"|) is a system of ordinary differential equations 
with constant coefficients. The tensors iYg, m and Gg^ m 
are given by one-dimensional integrals and are computed where 

I 



d t m x =-u) 1 m y 
d t m y =+LU 2 m x 



(28a) 
(28b) 



1 

w 2 



H cxt cos(9 - H ) + 2nM s (-N zz - N xx - (+N ZZ - N X3 >j cos(2#)) + 4^A/ s A 2 fc| m 

(29a) 

( -NZ (+Kz ~ N^x) cos(2#)) + 47rM s A 2 fc| m + 2nM s {l - 2G^J cos(20). 

(29b) 



This system of homogeneous equation has periodic so- 
lutions if the determinant of the characteristic system is 
equal to zero. This leads finally to the quadratic solution 

u)g=uiU2, (30) 

which defines the resonance condition for this geometry 



B. 3D micromagnetic simulation 

The dynamical susceptibility spectra of the small- 
est Py disk has been simulated by a micromagnetic 
code developed by S. Labbo^ and later by F. Boust 
and N. Vukadinovio^i. In this approach, the disk vol- 
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a) x"{x,y) at 8.2 GHz 
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FIG. 12: (Color online) Comparison between the 2D analyt- 
ical model (oblique lines) and the 3D micromagnetic simula- 
tion (horizontal spectra) for the prediction of the resonance 
fields of a D2 = 0.5 fj,m Py disk magnetized at 8h = 0°. The 
continuous lines are the locus of the resonance field of the 
radial modes (I = 0, m) predicted analytically (dashed lines 
correspond to the hidden modes (1 = 1, m)). 



ume is discretized by a regular cubic mesh of total size 
128 x 128 x 11, where each cube has an edge size of 3.9 
nm. The magnetization vector is assumed to be uniform 
inside each cell. This approximation is valid only because 
the cell size is smaller than the exchange length A = 5 nm 
for Py. Two 3D codes are used to calculate the dynamical 
response. For each value of the external field H cxt , the 
first code calculates the stable configuration of the mag- 
netization vector M(r) by solving the Landau-Lifshitz 
equation in the time domain. Hence, this code allows 
to incorporate in the simulation the small spatial depen- 
dence of the direction of M inside the sample volume in 
the quasi-saturated state. The second code computes the 
full dynamic susceptibility tensor \ from the lineariza- 
tion of the Landau-Lifshitz equation around the local 
equilibrium configuration. The used material parame- 
ters are identical to the ones measured in the extended 
thin film (7 = 1.849 10 7 rad-sec" 1 ^" 1 , 4ttM s = 9775 G, 
A = 10~ 6 erg/cm, a = 6 10" 3 ). 

At first, we shall compare the spectra calculated at 
Oh = 0° by the 3D micromagnetic simulation and the 
2D analytical model. Fig[l2] shows the amplitude of the 
imaginary part of the simulated susceptibility of the in- 
plane component x"xx at four different frequencies: 8.2, 
7.0, 5.6 and 4.2 GHz. A series of quantized modes is ob- 
served in the 8.2 GHz spectrum. The spatial distribution 
of the resonant modes in the mid-plane of the disks are 
shown in Fig ll3h for the three most intense peaks. The 
observed profiles correspond to the expected eigen-modes 
jf n {r), with t = and m = 0, 1 and 2 respectively. For 
comparison we have also plotted in Fig[l2] the locus of 
the resonances (see oblique lines) predicted by the 2D 
analytical model using the same parameters as in the 3D 



FIG. 13: (Color online) a) Images of the transverse suscep- 
tibility obtained through the 3D micromagnetic simulation 
for the first three modes at 8.2 GHz of a D2 = 0.5 /im Py 
disk magnetized at 6h = 0°. The thick lines indicate the 
nodes of precession, b) Image of the transverse susceptibility 
for the third mode (at H cxt = 8.8 kOe) of the 5.6 GHz spec- 
trum showing that the magnetization dynamics is not uniform 
across the thickness of the disk. 

model. The agreement between the two models is excel- 
lent for the lowest energy (highest field) modes. 

In FigfT3"b. the calculated transverse susceptibility of 
the third mode (at H cxt = 8.8 kOe) of the 5.6 GHz spec- 
trum is shown along the thickness of the disk. It shows 
that the magnetization dynamics is not uniform along 
the thickness, particularly in the center of the disk. This 
is due to the non-uniformity of the internal field along 
the thickness, as the disk is magnetized along its normal. 
This effect would get stronger if the thickness of the disk 
would be larger, and could eventually lead to an edge 
mode. In fact, the localization of the lowest energy mode 
at the top and bottom surfaces of Cu/Py/Cu sub- micron 
size disks, where the thickness of the Py layer was 100 
nm, was observed experimentally and calculated by full 
3D simulations of x"~- 

We have also performed micromagnetic simulations of 
the dynamic susceptibility at 6h 7^ 0°, using the same 
value for the parameters as above. FigQj] shows the sim- 
ulated spectrum at 8h = 5° for the four same excitation 
frequencies. The simulations are limited to the quasi- 
saturated field range when H cxt > 8 kOe. The highest 
frequency spectrum at 8.2 GHz shows two resonances. 
The spatial profile (not shown) of the lowest energy mode 
(located at 9.7 kOe on the 8.2 GHz spectrum) resembles 
J${r) described above. FigfTBTi shows the spatial profile 
for the second peak, at 9.2 kOe on the 8.2 GHz spectrum. 
The observed spatial profile breaks the axial symmetry 
due to the small tilt angle that the magnetization makes 
with the disk normal. This profile can be decomposed in 
the basis of the normal modes . The results is shown 
in Fig |15b . where we find that 

X"(r) « -^{^o°W-^i (r)-^ 2 (r-)| 

+ \{j2{r)+jKr)+Jt{r)}. (31) 
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0.9 1.0 

FIG. 14: (Color online) Comparison between the 2D analyt- 
ical model (oblique lines) and the 3D micromagnetic simula- 
tion (horizontal spectra) for the prediction of the resonance 
fields of a D2 = 0.5 fj,m Py disk magnetized at 6h = 5°. The 
continuous lines are the locus of the resonance field of the 
radial modes (I = 0, m) predicted analytically (dashed lines 
correspond to the hidden modes (1 = 1, m)). 



Hence, the normal modes in the presence of a small 
tilt angle Or are no longer the functions but rather 
a linear combination of them. We have also plotted in 
FigIH]thc 2D analytical prediction of the resonance fields 
as oblique lines. We find an important discrepancy be- 
tween the two models in the location of the resonance for 
the fundamental mode (almost 850 G apart), but also 
in the effective gyromagnetic ratio. In order to fit the 
3D simulated spectra with the 2D analytical model and 
9h as the adjustable parameter, a much larger value of 
8h ~ 10° has to be used in the analytical model than in 
the simulation. This shows the limits of the two main 
approximations of the 2D model: homogeneous dynam- 
ics along the disk thickness and unperturbed normal 
modes basis to compute the resonance fields. 



V. COMPARISON WITH THE 
EXPERIMENTAL DATA 

In this section, we want to analyze the three spectro- 
scopic informations that are i) the position ii) the ampli- 
tude and iii) the width of the resonance peaks measured 
experimentally on the Py disks of thickness 43.3 nm and 
of diameters D\ — 1.0 /im and Di = 0.5 /nm. 



A. Position 

As mentioned earlier, the values of the resonance fields 
are an important indicator of the effective internal field, 
which should not be disregarded by using the magnetiza- 
tion as a fitting parameter. We shall rely on the value of 




FIG. 15: (Color online) a) Image of the spatial distribution 
of the resonant profile for the 2nd peak at 9.2 kOe on the 
8.2 GHz simulated spectrum of Fig |14l The cartography is 
shown both in the axial and radial middle sections, b) The 
mode shown in a) is now decomposed in the Jm basis. The 
figure reveals here the spectral weights of the different eigen- 
modes. 
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FIG. 16: (Color online) Resonance field (blue dots) of the 
fundamental mode measured on the Di = 1.0 ^tm disk along 
with the analytically predicted angle dependence of H rea (u} a ) 
for different values of 8h (lines). Also indicated is the mea- 
sured resonance position in the extended film (red triangle) 
and the gyromagnetic ratio (dashed line): the shift with the 
8h = 0° line is due to finite size effects in the disk. 



M s measured by cavity-FMR on the extended thin film 
used for the fabrication of the disks and we shall assume 
that 4irM s = 9775 G is unchanged after the pattern- 
ing process. The cavity-FMR also provides the value of 
7 = 1.849 10 7 rad.sec^.G" 1 for our alloy composition. 
We first concentrate on the position of the fundamental 
mode (highest field mode) measured on the largest disk, 
of diameter D\ = 1.0 /mi (see Fig|TJ]). The blue dots 
shown in Fig |16l are the positions of the main peak put 
in a (uj s , -ff a ppi) diagram. Also on the figure (red triangle) 
is the resonance position of the uniform mode measured 
by the cavity-FMR, and the dashed line is the extrapo- 
lated locus of the mode as a function of frequency. We 
use the 2D analytical model presented in section HV Al to 
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fit the resonance fields. The results are shown in contin- 
uous lines. The two fitting parameters are the angle Oh 
and the offset field induced by the tip H a g. Changes of 
i?off corresponds simply to an over-all shift of the whole 
set of curves. At H Q s = Oe, the Oh = 0° line is al- 
ready shifted by almost 1 kOe lower in field (higher in 
energy) compared to the extended film. This is due to 
the cost in exchange and dipolar energies in the nanos- 
tructure. Increasing Oh in FigfTHl amplifies this shift but 
it also reduces the effective gyromagnetic ratio 7 G ff , which 
becomes frequency dependent. There is a unique value 
of 8h that reproduces the experimentally observed fre- 
quency dependence of j e g(H) and the best fit is obtained 
for Oh = 4.9°. Such a value for Oh is plausible, since our 
setup docs not incorporate a precise goniometer, and the 
angular orientation is difficult to tune on small size sam- 
ples. The fit also provides a value for H s = — 70 G. This 
value is close to the expected shift induced by the stray 
field of the tip at the distance s = 3.1 /mi of the disk. 

Using the fit values above, we have reported on FigfS] 
the locus of the higher order modes. The continuous lines 
indicate the cylindrically symmetric modes (I — 0) and 
the dashed lines indicates the modes with non-zero an- 
gular numbers i ^ 0, which should be hidden due to the 
homogeneous excitation field h. We find a good agree- 
ment with the data, except at the lowest frequency (4.2 
GHz) where a peak appears between the J® and J® 
modes. This new resonance is actually close to the lo- 
cation of the hidden mode. A possible explanation is 
that this resonance is a reminiscence of the broken sym- 
metry mode found in the simulation of the D2 = 0.5 
/xm disk at Oh = 5°. Such a mode, being a combina- 
tion of mainly J® and J® , would resonate between the 
m = and m = 1 mode. This mode becomes promi- 
nent at the lowest frequency because, as the excitation 
frequency decreases, the applied field becomes lower and 
the magnetic configuration becomes more sensitive to the 
in-plane component of the applied field, which favors the 
excitation of a mode having the symmetry of the in-plane 
component. 

The differences between the data and the 2D analytical 
model can be due to the approximations and assump- 
tions made in the model. One factor that affects the 
peak position is the value of the pinning condition. The 
expressions derived above in the 2D approximate model 
assume a total pinning of the SW at the periphery. As- 
suming no pinning at all at the periphery would move 
the peak position upward in field by about 50 G. It cor- 
responds to keeping the total pinning condition and in- 
creasing the disk diameter by some amount of the order 
of the film thickness^. Also, the radial component of the 
stray field of the probe is not taken into account in the 
model. This would accentuate the effects of Oh on the 
static configuration of M(r) in the sample. Other exper- 
imental uncertainties would be sufficient to account for 
the small discrepancy found in our analysis. We recall 
that the diameter of the disk used in the model is the 
one measured by SEM (see Fig[5]). A peripheral oxida- 



tion of the alloy can not be excluded, since there is no 
protective alumina on the periphery of the disks. Change 
of the disk diameter in the analytical model would shift 
the resonance peaks and affect the field separations be- 
tween them. Also, the separation between the probe and 
the sample, and the magnetic moment of the spherical 
probe, are known within 10%. This shows the limits of 
using the mechanical-FMR for a precise determination 
of the unperturbed FMR peak positions of a sub-micron 
size sample. 

Finally we have repeated this analysis on the smaller 
disk, with a diameter of D2 = 0.5 yum. Fitting the 
data with the analytical model gives 0h = 9.5° and 
H Q ff = —760 Oe. We have reported with straight lines 
on FigfTnithe results of the analytical model. The value 
of Oh found here is in disagreement with the previous 
finding, which is somewhat surprising since the two disks 
are located nearby, on the same substrate. Several ar- 
guments suggest that the true value of Oh is the value 
for the largest disk, i.e. « 5°. As shown in Fig[TJJ the 
analytical model tends to underestimate the effect of the 
angle as the disk shrinks in diameter, compared to the 
3D simulation. This affects the shift in the resonance 
field and the lowering of 7 ff, both underestimated by 
the 2D approximate model. A fit with this model will 
thus yield to a larger 0h value than the 3D simulation. 
If we now compare the experimental results of Fig llOl 
with the 3D simulation at Oh = 5° of FigfT4l we find 
that H g ss —1 kOe, which is a shift larger than the 
one expected, « —520 Oe from s = 1.0 fxm (see sec- 
tion III C 2[) . We emphasize here that the uncertainties 
mentioned above in the properties of the sample and the 
probe could translate in substantial errors in the 2D and 
3D analytical models. In fact, finite size effects in the 
D-2 = 0.5 /zm crucially depend on the exact properties of 
the sample. Moreover, due to the small separation be- 
tween the probe and the sample, the effect of the radial 
stray field from the probe is also more important than 
for the largest disk. 

B. Amplitude 

The amplitude of the modes indicates their coupling 
with the spectrometer. In the following we attempt to 
predict the relative amplitudes between the peaks. We 
will concentrate on the 5.6 GHz spectrum measured for 
the largest disk (cf. Figi]). 

Since the mechanical-FMR is using different schemes 
for the excitation and detection part, they need to be 
treated separately. We first compute the coupling to the 
excitation, which is identical to all FMR spectrometers. 
The microwave field being uniform at the scale of the 
sample, the coupling is simply given by the overlap inte- 
gral with the transverse susceptibility. In our notation, 
we find that the angle of precession is given by 

*■ = 3^ <32 » 
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FIG. 17: (Color online) Comparison between experiment and 
analytical model of the amplitude of the peaks for the disk of 
diameter 1.0 /im at 5.6 GHz. 

where we only consider the coupling to the (£ = 0, to) 
modes and where a is the damping coefficient. 

In our case, the force induced on the cantilever is given 
by 

/ AM z , m (r)3 zz (r, z + s + ^/2)d 2 rdz, 
Jv, 

(33) 

where the integral is the gradient of held along the 
z direction induced by the local variation of longitu- 
dinal component of magnetization inside the sample: 
AM, m (r) = \M s d 2 m J r l n (r) 2 /C^ rn and where 

S " (r ' Z) = ( r 2 +z 2)5/2 - ( r 2 + 2 2 )7 /2 ' ( 34 ) 

In Eq[33J we have used the fact that our probe has a 
spherical shape and can be viewed as a magnetic dipole 
TOsph placed at its center. 

We have reported in Fig[T7| a comparison between the 
measured spectra and the calculated peak shape using 
the analytical model and assuming that the line width is 
identical for all the modes. We find that the agreement 
with the data is excellent. In particular, we observe that 
the peak amplitude decreases by almost an order of mag- 
nitude between the J® and J® modes. This decrease is 
less pronounced in conventional FMR which measures the 
transverse component of the magnetization (see FigfT2|). 
This implies some caution when comparing the relative 
amplitude of the peaks found in the 3D simulation of 
the transverse susceptibility with the one found in the 
mechanical-FMR experiment. 

C. Line width 

Another important characterization concerns the 
width AH of the resonance line. The full circles and 



FIG. 18: (Color online) Frequency dependence of the line 
width measured on both the Di =1.0 disk (full circles) 
and the D2 = 0.5 fim disk (open boxes). The red triangle 
indicates the value measured by cavity-FMR on the extended 
film on Si. The mechanical-FMR data are renormalized by 
7 c jT. The dashed line corresponds to a damping coefficient 
a = 6 x 10" 3 . 

open boxes in FigdH] are the line widths measured as 
a function of frequency in the Py disks of diameters 
D\ = 1.0 fxm and D2 = 0.5 /im, respectively. All values 
have been renormalized by the effective gyromagnetic ra- 
tio 7 ff = du> s / dH &vp \. These line widths arc among the 
smallest reported for Py, confirming the excellent quality 
of the material used for this study. 

The resonance line width (half-width at mid-height of 
the absorption curve) of the Py polycrystalline magnetic 
layer is described^ as the sum of two contributions: 

AH = -4- + AH ICS (AH U , A9 H ). (35) 
27T1 

The first term relates to the intrinsic relaxation rate of 
the magnetization vector whereas the second term cor- 
responds to the inhomogencous broadening induced by 
a distribution of internal field, as well as by a distribu- 
tion of polarization angles, of characteristics widths AH U 
and A0h , respectively. AH U reflects the spatial inhomo- 
geneity of the field parameter H u whereas A9h refers to 
the distribution of the orientation Oh of the applied field 
with respect to the normal to the surface of the individual 
crystallites. 

The results obtained on the disks should be compared 
to the value observed by cavity-FMR on the extended 
thin film. In the cavity-FMR studies, the line widths 
measured on the film grown on the atomically flat Si sub- 
strate in the parallel and in the perpendicular geometries 
are equal (25 Oe) and reflect almost entirely the homo- 
geneous contribution. The measured value is represented 
by a red triangle in Fig(TS] It allows for an estimate of 
the intrinsic damping parameter (a = l/(2a;Ti)), found 
to be a — 7 10~ 3 . This value is an upper bound for a, 
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since the presence of a small inhomogeneous contribution 
in the line width would reduce the value of the homoge- 
neous part. The film grown on mica presents more inho- 
mogeneities than on Si, the results of which corresponds 
to a line width increase of 3 Oe. This is due to the even 
surface of the mica substrate, which bends over because 
of the strain between the different sheets of mica. 

From the frequency dependence of the line width of the 
large disk, we can estimate the damping coefficient in the 
nanostructure, found to be a = 6 10~ 3 (dashed line on 
Fig J18|) . in good agreement with the upper bound found 
on the thin film. We find that there is no inhomogeneous 
broadening in the nanostructure, which means that the 
small amount observed in the extended film deposited 
on mica is not relevant in structures confined at the sub- 
micron length scal o 10 ' 13 . ft also implies that there is no 
additional broadening induced by the mechanical-FMR. 
Finally, the frequency dependence of the line width at 
low frequencies and of the small disk is not linear. In 
fact, the magnetization configuration becomes not uni- 
form at low applied field, and the renormalization of the 
field line widths by 7 e g is not sufficient to recover the in- 
trinsic behavior. Moreover it was shown on the extended 
thin film that as the applied field is decreased, the in- 
crease of the angle between the equilibrium magnetiza- 
tion and the normal of the film increases the line width, 
through the inhomogeneous contribution associated 3 ^ to 
the distribution A8h in Eq|3"5l 



VI. CONCLUSION 

It was shown in this article how MRFM can be used 
to detect and quantitatively analyze the intrinsic FMR 
spectra of individual sub-micron size samples. To realize 
this objective, the coupling between the magnetic probe 
attached to the cantilever and the sample has to be op- 
timized. On one hand, it has to be as strong as possible 
to detect FMR in tiny samples. On the other hand, the 
inhomogeneity induced by the stray field of the probe 
has to be small compared to the internal dipolar field 
inhomogeneity in the sample, in order to detect the in- 
trinsic behavior of the latter. Using a magnetic sphere 
whose size is of the order of the disk-shaped samples en- 
ables to meet these two requirements. We demonstrate 
a 1000 spins sensitivity at room temperature. Using two 
approximate 2D and 3D models, we can understand the 
measured FMR-spectra on Py disks patterned out from 
the same extended film, whose characteristics are well 
known. It requires a good understanding of finite size ef- 
fects and of the homogeneous shift of the the SW modes 
resonance fields induced by the probe. 

Finally, we would like to summarize the main advan- 
tages of mechanical-FMR: i) its sensitivity to detect a 
single magnon excitation in a buried hybrid structure, 
e.g. below contact electrodes, ii) its versatility as a near 
field technique (i.e. only sensitive to the area directly 
underneath the probe), which allows spatial imaging of 



the magnetization dynamics, and iii) its ability to mea- 
sure the longitudinal component of the magnetization, a 
quantity directly linked to the damping. 

We arc greatly indebted to O. Achcr and A.-L. Adenot 
for their help and support. This research was partially 
supported by the ANR PNANO06-0235 and by the Eu- 
ropean Grant NMP-FP7 212257-2 MASTER. 



APPENDIX A: DEPOLARIZATION FACTORS OF 
A CYLINDER 

The analytical formula that have been used for the 
demagnetization tensor, TV, of a disk of radius R and 
thickness t are explicitly written in this appendix. The 
formula are actually derived from the published work of 
S. Tandoni^ but we have chosen to reprint them below 
because a couple of small typos remain in the original 
paper. The only assumption made here is that the mag- 
netization is homogeneous inside the cylindrical volume. 

Because of the axial symmetry, the values of the ten- 
sor TV arc better expressed in the cylindrical coordinates 
(r, z) and, using the notation of Ref4£, we introduce the 
reduced units £ = z/R,t = tj (2R), p = r/R: 

N zz (r, z )= + - {s C:r I {p, a_) + I (p, a+)} (Ala) 

N zr {r, z)=-\ {hip, <*-) - h(p, «+)} (Alb) 

N rr {r, z) = + - {s C:T I 2 {p, a_) + h(p, a+) - 2H r £h{p, 0)} 

- 4 {sc, T Io(p, a-) + I {p 7 a+) - 2H T . x I (p 7 0)} 

(Ale) 

where the notations a_ = |£ — r| and a+ = |£ + r| are 
respectively the distance (in reduced units) with the bot- 
tom and top surface of the cylinder. The function s and 
H design respectively the Sign and Heavisdc functions: 

1+1 if x < y 



H , 



— 1 else 



1 if x > y 
else 



The integrals U have the following expressions: 

kct 

I (p, a) =s ltP -A (/3, k) - - — -K(k) + H hp (A2a) 

Z ZTTy/p 

h(p,a)=-^-{(2-m)K(k)-2E(k)} (A2b) 
hM=^E(k)-(a* + p* + 2)^K(k) 



1 



si,p7r-oA (/3, k) 



(A2c) 
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where 



=fc 2 = 







arcsm 




(A3a) 
(A3b) 



K(k) and E(k) are the complete elliptic integrals of the 
first and second kind and Ao is the Heuman's Lambda 
function. 

The above expressions are valid everywhere in space. 
The magnetic field induction at every point in space (in- 
side or outside the sample's volume) simply obeys the 
formula: 

B z (r, z) = H CKt + 4ttM s {6(r, z) - N zz (r, z)} , (A4) 

where O is a function equals to 1 inside the cylindrical 
volume and otherwise: 



9 around the y direction. If 9r is the angle that the ex- 
ternal magnetic field makes with the normal of the disk, 
then 9 is implicitly defined by the equilibrium condition: 



iJext sin(0 -9 H ) + 2ttM s N xx sin 29 - 2irM s N zz sin 29 = 

(Bl) 

In all the expressions above, the tensors are expressed 
in the Cartesian frame of the disk (x,y, z). This applies 
for the demagnetizing factors N in appendix fXl but also 
for the matrices: 





1 



1 





(B2) 



(B3) 



Q(r,z) 3 



1 if r < R and \z\ < t/2 
else. 



APPENDIX B: LINEARIZATION IN THE LOCAL 
FRAME 

The notations have been defined in Figllll (x, y, z) 
is the Cartesian frame along the principal axis of the 
disk, with z oriented along the normal and (£, y,C) the 
Cartesian frame of the magnetization dynamics, with £ 
along the equilibrium (or effective magnetic field) direc- 
tion. Both frames are related by a rotation of an angle 



Solving Eq. (|24p requires to write down the different 
tensors T(= N, G) in the local frame of the magne- 
tization. This is achieved through the transformation 
RgTRg, where 01$ is the rotation matrix between the 
(x,y,z) and (£,y,C): 



(B4) 



cos 9 0— sin ( 

1 
sin 9 cos 9 
and the subscript t stands for the transpose. 

This allows us to find the expression Eqs. (f2"9")) . where 
we have used symmetry arguments to impose N xy = 
A^I=0. 
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